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Abstract 

In this work the most spectacular phenomena occurring during the formation of 
a star are briefly reviewed: accretion through a rotating disc of matter and outflow 
through the poles of the new stellar object. Magnetic fields have been proposed 
to be principally responsible for the coexistence of these opposed mechanisms of 
accretion and outflowing. According to different models, the magnetic fields are 
either twisted or stretched by the accretion disc, allowing the formation of polar 
channels through which part of the accreting matter can escape. The high degree of 
coupling between the physical parameters describing the young stellar object (e.g. 
circular velocity of the disc and the central object, viscosity, strength and freezing 
of magnetic fields, etc.) makes detailed understanding of the interaction between 
accretion discs and outflows very difficult. The models presented here provide only 
a partial answer to this difficult problem. 



1 Introduction 



Observing a main sequence star, like our Sun, but much further away with an optical 
telescope will not impress us much, since it will appear as an unresolved point source of 
light, without any particular feature. Main sequence stars are in fact objects in the most 
stable and quiet stage of their life. A slightly higher degree of activity (e.g. flares) can 
be observed at other wavelengths (e.g. radio or X-rays): synchrotron emission from the 
very hot plasma in the corona of the star is one of the most evident hints for the existence 
of strong magnetic fields rooted in the star. Nevertheless this activity is not comparable 
to the one through which the star had to go before reaching the quiet days of the main 
sequence. First the collapse of a molecular cloud, which can be of several light years in 
size. Then the accretion of matter through a disc accompanied by outflows (see Fig. |l]). 
Finally the star reaches the stage of the ignition of nuclear fusion in its core which leads 
to emission at optical wavelengths, which will characterize its life on the main sequence. 

I concentrate on the second phase mentioned above and explain, on the basis of both 
on theory and observations, the coexistence of two apparently opposite phenomena as 



*Astrophysical Dynamics 1999/2000, Alessandro B. Romeo (Ed.), Onsala Space Observatory, 2000. 
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Figure 1: Single wavelength Hubble Space Telescope (HST) images do not give infor- 
mation about the dynamics of young stellar objects but they are able to constrain their 
geometry. Here is shown HH 30, a T Tauri star, with a disc and two opposite polar 
outflows. Picture taken from the NASA archive. 



accretion and outflow of matter. To do this I will first introduce briefiy (in section 2) 
the physical principles which are needed for such a study. I will then review a sample of 
theoretical models and the results of recent simulations (section 3). In section 4 I will 
draw the best links between the model predictions and the observed phenomena. As a 
summary, in section 5 I will try to delineate the further work in the field. The last section 
summarizes the discussion occurred at the end of the oral presentation. 



2 Basic Equations 

Magnetohydrodynamics is the study of the motions of a fluid in the presence of magnetic 
fields. We consider the material accreting and outflowing from a star as a fluid because 
the mean free path / of a particle is much smaller than the length of interest L (in our case 
/ ~ 10^^ cm and L ~ 10^^ cm. We therefore will need a set of equations which include 
both the dynamics of a fluid and the presence of magnetic flelds. 

First the conservation of the mass from the hydrodynamic equations: 

^p + V(pv)=0 (1) 

We then need an equation for the conservation of momentum: 

P[^v + (v • V)v] + Vp - pg = ^ (V X B) X B (2) 

Notice that the term on the right marks the presence of the magnetic fields with the 
magnetic force, p represents the pressure and g the gravitational potential (which can 
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have whatever shape we wish). 

The energy equation looks hke the following: 



+ V ■ (Ev) = -gv + i-r/ I V X B |2 (3) 

where rj corresponds to the electric diffusivity. 
Finally we need field equations: 

d 

— B - V X (v X B) = -V X (?7V X B) (4) 
V-B = (5) 

This set of equations needs a few lines of explanation. First, in magnetohydrodynamics 
(MHD) we assume local neutrality. Because of the very low density, the charged particles 
of both polarities (electrons and ions) are supposed to fill the same space equally. This 
does not rule out the existence of different drift velocities between the two species. It 
means that the currents exist but the electric fields can be neglected. Hence the Lorentz 
force looks like the following: 

-J X B = —(V X B) X B (6) 

Since the electric fields are considered small enough to be neglected, their changes in time 
are also negligible. So we have 

4:71 Id 

where the last term will be equal to 0. In addition, if we rule out all dissipation effects 
[rj = 0) because of the very low density of the plasma, we have the right hand side of 
equation equal to 0. In that case we talk about field freezing, which means that the 
field lines follow every movement of the matter in which they are embedded, producing 
changes in the geometry of the magnetic field and in the strength of the magnetic flux. 

The solution of this system of equations is a great challenge in both theoretical mod- 
elling and the numerical computation. It will trace the magnetic field lines around young 
stellar objects and meet the observations which have been carried out with all kind of 
techniques and at all wavelengths. In the next section we will have a look at the models 
which give the best fit to the observed data. 



3 The Models 

From the literature it is possible to recognize two major schools in modelling the magnetic 
fields around young stellar objects. The first (e.g. Shu et al., 1988) considers the outfiows 
as centrifugally driven winds; the second stresses the twisting of the magnetic field lines. 

3.1 Magnetocentrifugally Driven Outflows 

Shu et al. (1988, 1991, 1994) proposed a model where the outfiows detectable around 
young stellar objects are driven basically by the centrifugal force produced by the fast 
rotating central object. The latter is supposed to rotate nearly at break up speed. The 
poloidal magnetic field lines are squeezed towards the surface of the star, producing big 
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Figure 2: The model of magnetocentrifugally driven outflows by Shu et al. (1988). Note 
that, as a particle recedes from the star along a field line, the gravitational force 
becomes weaker, whereas the centrifugal force Fc becomes stronger. 



loops expanding along polar directions. The accretion disc drives the matter towards 
the star. Some of the matter travelling onto the disc surface to the star can leave the 
plane of the disc and fly along the squeezed field lines. In Fig. || it is possible to see the 
geometry of the field lines close to the central object and how the centrifugal force Fc 
wins over the gravitational attraction F^ the further the particle travels along the field 
line (Pelletier et al., 1992). The problem with this model is the collimation of the outflow. 
Observations show in some cases evidence for a very efficient mechanism of collimation of 
the matter flying away from the central object. This effect is not reproducible with this 
model. An example of an application of the model is given in Shu et al., 1994, on SU 
Aur a classical T Tauri star. = 2.25Mq, = 3.6Rq give a mass flow from the disc 
Md = 6 X 10~^ Mq/uv and a magnetic field strength B^, = 300G. The fiow velocity will 
be around 200 km/s. 

3.2 Twisting the Magnetic Fields 

In the following models the mechanism that drives the outfiows is not a wind created by 
the centrifugal force. The forces acting as launching the matter away from the central 
object are considered to originate only by the twisting of the magnetic field lines, which 
is able to open polar channels and draw matter from the accreting disc towards the 
interstellar space. 

The first of these was proposed by Uchida & Shibata, 1985. It considers an interstellar 
magnetic field, frozen into the interstellar medium. Due to the collapse of the original 
molecular cloud the magnetic field lines are constrained to come closer to the central 
object, and because of the rotation of the accretion disc they will be twisted. The result 
is shown in Fig. ^ The twisting of the field leads to the appearance of a J x B force: J is 
radial due to infalling material, B leaves its initial configuration parallel to the rotation 
axis and, because of the twist due to the rotation of the disc, lies more and more in the 
plane of the disc. The J x B force is therefore oriented perpendicular to the disc and is 
able to produce large scale weakly collimated outfiows. The collimated hot optical jets. 
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Figure 3: The model by Uchida & Shibata (1985). The three panels show a progressive 
zoom into the star formation region. Note in the right panel the closed and opened field 
lines. 



on the other hand, are considered to have their origin closer to the star: material falling 
to the surface of the star can be blown off in the polar direction along the "open" field 
lines (these are not strictly open, they in fact reconnect with the interstellar field). 

The second model comes from Newman et al. (1992) (Fig. ^). The important assump- 
tion here is J x B = 0, which indicates a force free plasma. The Poynting theorem applied 
to this condition shows that because of the differential rotation of the disc (footpoints of 
the poloidal field travel apart form each other) the energy of the field has to increase. 
This energy enhancement causes the field lines to "infiate", especially in the direction 
of the poles. Material can then escape along the infiated field lines, away from the disc 




Figure 4: The Newman et al. (1992) model shows how differential rotation can twist the 
magnetic field lines, forcing them to expand in the polar direction. The time sequence 
starts from the left drawing. 
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surface. Simulations by Lovelace et al. (1995) show that very little time is needed in 
order to reach a reasonable field infiating, which means that the phenomenon of outfiow 
is present at the very beginning of the formation of the star, as soon as the accretion disc 
is present. 

The third model was proposed by Henriksen & Valls-Gabaud (1994) and is called cored 
apple. It is explicitly said in the article that this model cannot reproduce any observation 
with high fidelity. Nevertheless it serves to show that no possibility can be ruled out. In 
particular their simulation shows, beside a known equatorial- infall-/polar-outfiow-model, 
the opposite scenario of having a polar infiow (better called infall) and an outfiow in 
the equatorial plane as (see Fig. ||). This mechanism is supposed to appear at the very 
beginning of the star formation, i.e. during the collapse of the molecular cloud, when yet 
no disc-like geometry has been created from the accreting material. Only a thick torus is 
present. The engine of the outfiow is only the rebounding of the free falling material. It is 
even possible for some of the gas to miss the central object and to be directly conducted 
outwards. The model is however not able to reproduce field energies, pressures and fiow 
velocities measured by observation, especially when the distance to the central object 
becomes greater than 100 AU (tests of the model have been carried out on VLA 1623, 
see Henriksen & Valls-Gabaud, 1994). 

There have also been several attempts of modelling time-dependent accretion and 
outfiow, see e.g. Goodson et al., 1997, 1999a, 1999b. They basically assume a radially 
oscillating disc, which periodically squeezes the poloidal magnetic field (field configuration 
of Newman et al., 1992), which infiates and allows material to fiow out in the form of 
"bullets". The model is in good agreement with observations at different wavelengths 
of HH 30 (see Fig. |I|) and DG Tau, two classical T Tauri stars. Further attempts of 
describing the material ejection time dependence claim to trace both the phenomena 
around a black hole and around a protostar (Ouyed et al., 1997). 

The latest models are able to combine the existence of a fast, hot and narrow outfiow 
(near to the rotation axis) with a slow outfiow, colder and closer to the disc. The first 
would be driven by the dynamo action of the disc which gives raise to an azimuthal field, 
the second would be centrifugally driven. The former would be responsible for carrying 
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Figure 6: The emission spectrum of the CO transition J = 2 — 1 towards the outflow in 
NGC 2264G. The upper panel shows the emission of the redshifted lobe of the outflow, 
the lower panel the emission of the blueshifted (Lada & Fich 1996). 



material away from the disc, the latter to transport angular momentum outwards, along 
the disc (Brandenburg et al., 2000). 

4 Observations: The Reality 

From the pictures of the HST we get a static image of a forming star (Fig. |l]). All 
information of the dynamics in that region is completely missing. In order to supply 
this information it is important to combine different techniques. In particular, I will 
present some results obtained from spectrometric studies of those regions at millimeter 
wavelengths. The reason for using radio observations instead of optical ones is simple: 
optical radiation is absorbed by the great amount of dust present around the central 
object; radio waves on the contrary get through and carry very important information 
about the interior of the cloud. Infrared observations are also important, since the dust 
scatters the [/^-light from the star to longer wavelengths producing a red excess, i.e. a 
higher flux in the infrared than would be expected from the star. This excess radiation 
gives us information about the amount of dust around the forming star. 

An observed spectral line will often be shifted in frequency, because of the relative 
radial velocity between us and the observed object (Doppler-effect). From the amplitude 
of that shift and from the shape of the spectral line it is possible to draw conclusions 
about the relative motion of the matter, e.g. if it is moving towards us or away from us. 
In Fig. ^ a typical spectrum is shown. Using an interferometer the data can be mapped in 
both space and velocity to study the source dynamics. Clearly noticeable is the velocity 
distribution of the matter at different distances from the central object (Fig. |^). 

But even more useful for the study of the magnetic fields is the measure of the po- 
larization of the light. By observing at millimeter and submillimeter wavelengths it is 
possible to measure the amount and direction of polarization as e.g. shown in the paper 
by Tamura et al. (1999). Dust particles aligned with the magnetic fields scatter UV- 
light from the star producing a slight polarization. In Fig. § one can see the orientation 
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Figure 7: From the spectrum above it is possible to map the emission spatially: every 
slide corresponds to a range relative velocities, showing how the material is distributed 
both in space and velocity (Lada & Fich 1996). 



of the polarization (thick line) superimposed on CO contours which traces the outflows. 
Realizing that the magnetic field is oriented along the disc, one can conclude that the 
toroidal component of the magnetic field is dominating the configuration around the two 
objects. It is also conceivable that determining the age of young stellar objects will be 
possible through measuring of the polarization: the more alignment with the accretion 
disc is detected, the older is the YSO (as time goes by the configuration is geometrically 
more defined and the magnetic field more twisted in the toroidal direction). A discussion 
about the possible correlation between "order" in the magnetic field and age is presented 
by Greaves et al. (1997). 

The choice between one theoretical model and another is not always clear. Fig. ^ 
shows a plot of polarization against viewing angle. The latter is calculated based on 
geometrical considerations (see Greaves et al., 1997). The clear correlation between the 
polarization angle and the angle of view is striking and this suggested to the authors dif- 
ferent conclusions about the geometry of the star- disc-outflows system. From the diagram 
in Fig. ^ one can see that when the outflow lies in the plane of the sky (B = 90°) the 
magnetic field lies perpendicular to the outflow direction. The opposite happens when 
the outflow is oriented close to the line of sight: the magnetic field seems to be parallel 
to the outflow. The authors test against the last three models of section 0. Though a 
definitive choice between models is not possible they can nevertheless possible to rule out 
one of them in some cases. 
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Figure 8: Two examples of a measurement of the polarisation of the mm or sub-mm 
continuum emission coming from young stellar objects (in this case two classical T Tauri 
stars). The thick hne indicates the orientation of the magnetic field, which lies perpendic- 
ular to the outflows shown in the CO. One can conclude that the toroidal component of 
the magnetic field dominates the configuration around the YSO, indicating an advanced 
state of evolution (the magnetic field lines have been twisted for long time) (Tamura et 
al. 1999). 
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Figure 9: This diagram shows the correlation between A© (angle between the magnetic 
field and the outflow direction) and 6;os (angle between the outflow and the line of sight 
to the observer, estimated geometrically. From this diagram it is possible to determine 
the geometry of the object (Greaves et al. 1997). 



5 What Remains to Be Done 

In my opinion, the greatest challenge of the next generation of instruments (ALMA for 
the submillimeter, optical and infrared interfcromctry, etc.) is to achieve higher spatial 
resolution in order to be able to zoom into the deepest regions of the formation of a star. 
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Polarization measurements are also well suited to study the presence and the evolution 
of magnetic fields in star formation regions. All this will be a great contribution to the 
modelling of both jets from black holes and outflows from YSO's. In active stars (pre main 
sequence), the tracing of the magnetic fields can also be done by observing synchrotron 
emission from the plasma of the corona. These studies give insight into the existence and 
evolution of the corona in early phases of stellar evolution, in order to better understand 
the origin of e.g. the solar corona. An ongoing project tries in fact to detect T Tauri 
stars at very high resolution through synchrotron emission, with the aim of tracing the 
distribution of the hot plasma and consequently of the magnetic field. 

6 Questions and Discussion 

• Student: Where does the outfiowing matter go? Does it come back on the star 
later on? 

M. Pestalozzi: The matter is blown away from the star and will in principle not 
participate in the formation of the star where it originated. Eventually it might 
be part of the formation of another star in the neighbourhood. When stars form 
in clusters the star formation rate is triggered exactly by this kind of phenomenon, 
which spreads out material into the interstellar space. 

• A. Romeo: How does accretion occur? 

M. Pestalozzi: A single particle moving in a central field has a constant angular 

momentum. This fixes the characteristics of the orbits that the particle can have 
(the shorter the radius, the higher the rotation velocity). In order to fall towards 
the central object the particle has to lose angular momentum and energy, which in 
a YSO occur through the frequent collisions of gas particles. Research is in progress 
to understand how these dissipative phenomena are correlated. 

• A. Tappe: Jets from a black hole and outflows from a YSO: are they the same 
phenomena? Will it be possible to join them in the same class? 

M. Pestalozzi: In the simulations by Ouyed, Pudritz & Stone (1997) the prop- 
agation of jets and outflows does not differ, except for the scale. However YSO 
and black holes are physically very different objects, especially close to the center. 
The uniflcation of the theories describing outflows or jets is probably not possible 
at scales comparable to the size of the central object. Further out the morphology 
of the dynamics of the outflowing material have actually similar characteristics in 
both type of objects. 

• M. Thoiricisson: How does polarization of light occur? 

M. Pestalozzi: Light coming from the star scatters on the dust particles which lie 
around the star (either in the disc or in general in the cocoon of the forming star). 
Magnetic fields on the other hand have the capacity of orienting the dust particles 
along the field lines, according to their magnetic properties: every elongated particle 
will feel a magnetic momentum, which will force the particle to assume the less 
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expensive position in terms of energy. In this situation, light which scatters with 
oriented dust will orient its fields, producing a polarization along the major axis of 
the dust particles. 

• M. Thomasson: Where is the force that drives the matter out of the plane (e.g. 
in Newman et al., 1992)? 

M. Pestalozzi: In the model by Uchida & Shibata, 1985 it is the J x B force which 
increase with the twisting of the field lines by the rotation of the disc. In Newman et 
al., 1992, material follows the field fines which have been infiated by the stretching 
of the field fines due to differential rotation of the disc. 

• A. Tappe: Open field lines: is this possible? 

M. Pestalozzi: What is meant by open field lines are not lines which will never 
close, opening the possibility to the existence of magnetic monopoles. These lines 
are supposed to reconnect with the fines of the interstellar magnetic field. They are 
virtually open for the central object we are considering. 

• A. Tappe: Is it possible to tell which model is right and which is wrong? And if 
yes, how? 

M. Pestalozzi: No. There are better models and worse ones. The whole thing is 
about how well a certain model fits the observations. Some of the presented models 
are better for the explanation of some observations but not of others (see Greaves 
et al., 1997). 

• V. Minier: You are speaking about low-mass star formation. What about high- 
mass star formation? Is any model possible for it? 

M. Pestalozzi: I was actually expecting that question. In my opinion that I have 
been cultivating during the preparation of this seminar, there should not be any 
difference between the formation of a low-mass star and a high-mass star, except 
for the scale of the phenomena. The problem comes from the time scale and the 
chance to observe a star in this phase: high-mass stars are supposed to form and 
evolve quicker than low-mass stars. This already puts a constraint on the chance 
of finding a high-mass star exactly in the evolutionary phase we would like. In 
addition, there are many more low-mass stars than high-mass ones. In any case, 
basing our consideration on the existent literature, I suggest to look at a high mass 
star formation applying Henriksen's model. Since it seems to be difficult for a 50 
Mq star to keep intact discs or other nice features, the model of outflows created 
only by the accretion seems to me the most appropriate to describe the early phase 
of evolution. 

• A. Tappe: The collimation is proportional to the twisting, or better to the velocity 
difference between the rings where the footpoints of the magnetic field lie. Is that 
correct? 
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M. Pestalozzi: Yes. It is also well presented in the article to which I was referring. 
Actually in the proposed model the authors needed to turn the simulation only 
about a bit more than one turn in order to get those very pinched field lines. 

• D. de Mello: On the HST viewgraphs you showed at the beginning a lengthscale 
is indicated. Could you compare it to the one you get observing with VLBI? 

M. Pestalozzi: In fact with VLBI observation we are able to zoom in by about 
two orders of magnitude inside the objects we see in the HST picture. This seems to 
do all simulations about the magnetic fields useless and a comparison between HST 
images and VLBI data nonsense. In fact the simulations try to reveal real mechanism 
that starts the outfiows, almost at the stellar surface. What we see with the HST 
images is a constraint on the geometry at a larger scale. The models must then be 
able e.g. to coUimate the outfiows within the distance form the central object that 
HST images indicate. Also the accretion mechanisms must obey those geometrical 
boundary conditions. In such complex problems of MHD in YSOs where all effects 
seem to couple (viscosity, accretion, outflow, magnetic flelds, rotation, diffusion, . . . ) 
it is important to have at least some kind of boundary conditions on the geometry 
on which you can build a model. 
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